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Abstract. Enhanced production of events with almost equal-sized fragments is experimentally revealed 
by charge correlations in the multifragmentation of a finite nuclear system selected in 129 Xe central colli- 
sions on nat Sn. The evolution of their weight with the incident energy: 32, 39, 45, 50 AMeV, is measured. 
Dynamical stochastic mean field simulations performed at 32 AMeV, in which spinodal instabilities are 
responsible for multifragmentation, exhibit a similar enhancement of this kind of events. The above ex- 
perimental observation evidences the spinodal decomposition of hot finite nuclear matter as the origin of 
multifragmentation in the Fermi energy regime. 

PACS. 2 5.70.Pq - 2 4.60.Ky 



1 Introduction 

When enough energy is brought into a nucleus, it breaks 
up into smaller pieces: this is called multifragmentation 1 , 
The knowledge of this process, experimentally observed 
for many years, was recently strongly improved with the 
analysis of experiments performed with powerful detection 
arrays [5||l| | 5||^| |7llc1l51llUlllllll2| . As a theoretical frame- 
work for describing this phenomenon, the analysis of the 
bulk dynamics of nuclear matter, under various condi- 
tions of density and internal energy, has found arguments 
that multifragmentation occurs when nuclear matter has 
expanded through the spinodal region of negative com- 
pressibility |13|. In this area of mechanical instabilities, 
covering a large part of a liquid-gas type coexistence do- 
main, the irreversible growth of local density fluctuations 
leads the system to the separation into two phases: the 
spinodal decomposition. For the liquid part primary frag- 
ments of nearly equal sizes should be favoured, in relation 
to the wavelengths of the most unstable modes (T3|. Ef- 



a Present address: Cyclotron Institute, Texas A&M Univer- 
sity, College station, Texas 77845, USA 

Permanent address: Institut de Physique Nucleaire, 
IN2P3-CNRS, F-91406 Orsay cedex, France. 
Correspondence to: borderie@ipno.in2p3.fr 



fects like beating of different modes, coalescence of nascent 
fragments, secondary decay of the excited fragments and, 
above all, finite size effects are expected to deeply blur this 
simple picture [15116) . However, evidence for an enhance- 
ment of events with nearly equal-sized fragments was for 
the first time observed in 32 AMeV 129 Xe + nat Sn central 
collisions leading to a fused system which undergoes multi- 
fragmentation, by using a charge correlation method |17| . 
By means of this model independent method, recently im- 
proved |18j , we are investigating here the evolution of the 
equal-sized fragment partitions with the increase of the 
incident energy up to 50 AMeV, in selected samples of 
experimental fused events concerning the same system. In 
the same framework, for collisions at the lowest energy, the 
predictions of the 3D stochastic mean-field simulations are 
successfully compared to the experimental results; these 
simulations take into account the dynamics of the most 
unstable modes in the spinodal region. 

The paper is organized as follows. In section 2 we firstly 
present the experimental set-up, including the detector ar- 
ray specifications and operating conditions during exper- 
iment. Secondly, we recall the criteria allowing to select 
the experimental events corresponding to fused systems 
in central collisions. The general features of the related 
experimental data are evidenced in section 3. In section 
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4, the charge correlation method with some improvements 
and developments is described; results of stochastic mean 
field simulations are used to exemplify and to compare 
the potentiality of the up-dated version to the original 
method. The experimental results, obtained at different 
incident energies, are then shown; they are discussed and 
interpreted in section 5. Conclusions of this study are 
drawn in section 6. 



2 Experimental selection of a finite piece of 
nuclear matter 

Heavy ion collisions at relative velocities comparable to 
those of the nucleons in the nucleus may provide test 
pieces of nuclear matter at moderate temperatures which 
are good candidates to undergo a nuclear liquid-gas type 
phase transition. Indeed, "fused" systems reaching sizes 
up to a few hundreds of nucleons may be obtained in the 
laboratory, by colliding accelerated heavy projectiles with 
heavy targets. By means of 4w detection arrays of high 
granularity, it is possible to completely study their dis- 
assembly. Such events are expected to appear with high 
probability for central collisions, but one major problem 
is to select them among reactions with a dominant binary 
character |19j . 

2.1 Experimental procedure 

The 129 Xe + ™*Sn system was studied with the 4ir mul- 
tidetector INDRA, operating at the GANIL accelerator. 
A thin target of natural tin (350 fig/cm 2 ) was bombarded 
by 129 Xe projectiles at five incident energies: 25, 32, 39, 
45 and 50 AMeV. Low intensity beams (~ 3 x 10 7 pps) 
were used to keep the random coincidence rate below 10 . 
The low target thickness allows slow fragments to escape 
the target. A trigger based on the multiplicity was cho- 
sen, requiring at least four modules firing. In the off-line 
analysis, events having a multiplicity of correctly identi- 
fied charged particles inferior to the experimental trigger 
condition were rejected for reasons of coherency. The most 
peripheral collisions are thus eliminated. 

INDRA, which is described in detail in (20II21| . can be 
viewed as an ensemble of 336 telescopes covering about 
90% of the 47T solid angle. The detection cells are dis- 
tributed amongst 17 rings centred on the beam axis. Low- 
energy identification thresholds (from rj 0.7AMeV for Z = 

3 to w 1.7AMeV for Z = 50) and large energy ranges were 
obtained through the design of three layer telescopes, com- 
posed of an axial-field ionization chamber operated at 30 
mbar of C^Fg, a 300 //m silicon detector and a CsI(Tl) 
scintillator, thick enough to stop all emitted particles, cou- 
pled to a phototube. Such a telescope can detect and iden- 
tify from protons between 1 and 200 MeV to uranium ions 
of 4 GeV. Past 45°, where fast projectile-like fragments 
are no longer expected, the telescopes comprise only two 
stages, the ionization chamber operated at 20 mbar and 
the scintillator. Finally the very forward angles (2 - 3°) 
are occupied by NE102 - NE115 phoswiches. 



A charge resolution of one unit was obtained for the 
whole range of atomic number of fragments identified through 
AE - E method in the Si - CsI(Tl) couple. For CsI(Tl) 
scintillators a better understanding of the light response 
was obtained and the contribution of 6 — rays generated 
by the incoming heavy ion was taken into account |22|. 
A direct consequence was the identification of fragment 
with a resolution of one charge unit up to Z = 20, and a 
few charge units for the heaviest fragments, in ionization 
chamber - CsI(Tl) telescopes [SB]. The exact identification 
of fragments up to at least Z = 20 reveals essential for the 
charge correlation studies which are the aim of this paper. 

2.2 Selection of single source events 

A two step procedure has been used to isolate fused sys- 
tems. The first step was to keep the events for which 
a quasi-complete detection of the reaction products has 
been achieved. Significant fractions: > 77% of the charge 
of the system, Z sys = Z pro j + Z tar g, and > 75% of the 
beam momentum, P pro j , in the exit channel were required 
to be measured for every event. INDRA does not permit 
isotopic identification for fragments and does not detect 
the neutrons. For this reason, the momentum used here is 
calculated from the product of atomic number Z and ve- 
locity component in the beam direction v z : P to t = SZv z 
and normalized to the incident (projectile) momentum: 
ZprojVproj. In the second step we used the flow angle 
{Ofiow) selection [24125119] . This global variable is de- 
fined as the angle between the beam axis and the pre- 
ferred direction of emitted matter in each event. It is de- 
termined by the energy tensor calculated from fragment 
(Z > 5) momenta in the reaction centre of mass. Fused 
events have no memory of the entrance channel and should 
be isotropic while binary dissipative collisions arc focused 
at small 0fi ow - Thus, by selecting only large flow angles, 
fused events can be well isolated. In ^U] the minimum flow 
angle chosen for 32 AMeV collisions was 70°. It was 60° 
for 50 AMeV collisions [20] • In this paper we chose 60° 
for all energies, to get enough statistics (at least 30000 
events) without degrading the properties of a compact 
fused source. 

The present selection corresponds to measured cross 
sections decreasing from 56 to 19 mb when the incident 
energy goes from 25 to 50 AMeV (see table QJ. By tak- 
ing into account detection efficiency and biases due to the 
selection (quasi-complete events and flow angle selection) 
the total cross section for the formation of compact fused 
systems is estimated to decrease from 250 to 85 mb be- 
tween 25 and 50 AMeV [27]. 

3 Global properties of selected fused events 

Before presenting and discussing charge correlations, which 
deal with information inside events, it is useful to have an 
overview of the related inclusive properties of the selected 
events, like the multiplicity and charge distributions. En- 
trance channel conditions for the Xe+Sn system of total 
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Table 1. Exit channel average values of the total and fragment (Z > 5) multiplicity, total charge emitted in fragments and 
fragment atomic number, in events from 129 Xe + nat Sn central collisions, as functions of the entrance channel characteristics: 
incident and available energy. The numbers in parentheses are the standard deviations of the distributions. The measured cross 
sections are given in the last column. 



E 




(M tot ) 


(M f ) 


(Abound) 


(Zf) 


&meas. 


AMeV 


AMeV 










mb 


25 


6.24 


18.9 (2.7) 


3.56 (1.04) 


67.3 (6.8) 


17.6 


56. 


± 8. 


32 


7.99 


23.8 (3.0) 


4.13 (1.17) 


57.6 (7.0) 


13.4 


26. 


± 4. 


39 


9.73 


28.5 (3.1) 


4.41 (1.20) 


49.6 (7.4) 


11.0 


18. 


± 3. 


45 


11.23 


31.5 (3.3) 


4.42 (1.19) 


44.0 (7.5) 


9.7 


19. 


± 5. 


50 


12.48 


34.2 (3.3) 


4.31 (1.18) 


39.7 (7.6) 


9.0 


19. 


± 3. 



mass A sys = 248 and charge Z sys — 104 are summarized 
in Table □ 

3.1 Multiplicity and bound charge 

Fig. ^ shows multiplicity and total charge (Abound = 
12z>z distributions normalized to the number of 

events, at the five incident energies. The total multiplicity 
distributions, when all the charged reaction products are 
considered, Z > 1, are shown in fig.^i. They all display a 
gaussian shape; the most probable multiplicity values and 
the standard deviations increase with the available energy, 
as shown in Tabled At 50 AMeV the average multiplicity 
reaches one third of the total charge. The distributions of 
the sum of the charges of all the reaction products, plot- 
ted in fig.n?), for different energies, are practically identi- 
cal, nicely testifying about the negligible variation of the 
detection efficiency of INDRA when the incident energy 
significantly varies. The common limitation at 80 charge 
units is demanded by software as a criterion for the com- 
pleteness of event detection. These two plots, considered 
together, speak about the increasing number of reaction 
products from central collisions with the incident energy. 

The evolution of the reaction product multiplicity dis- 
tribution, at all energies, when another lower charge limit 
is imposed, is presented afterwards in the left column of 
the same figure. As compared to fig.^. (Z > 1), the split- 
ting of the energy dependent distributions is diminishing 
if Z > 3 (fig. ; the curves become completely superim- 
posable, except at 25 AMeV, for Z > 5 (fig. [Q) and split 
again, but in the opposite sense, if Z > 8 (fig.^i). Deexci- 
tation simulations for primary fragments having the same 
N/Z ratio as the total system, performed by means of the 
statistical code SIMON [2Hj at 3 AMeV excitation ener- 
gies, have indicated small evaporation rates of Li and Be 
isotopes, but negligible rates of heavier nuclei |2H]. This 
average value of 3 AMeV for the fragment excitation en- 
ergy was deduced from fragment-particle correlations for 
the same system between 32 and 50 AMeV incident en- 
ergy [3HII2T] . Moreover the daughter nucleus atomic num- 
ber is, on average, one charge unit smaller than that of 
the parent fragment, with a standard deviation of about 
one unit. It may be thus inferred that the distribution of 
the number of primary fragments with Z > 5 is not dra- 
matically modified by the secondary decay. More precisely, 



as indicated in Table the average number of fragments 
with Z > 5 reaches a maximum for an incident energy as 
low as 39 AMeV and starts decreasing above 45 AMeV. 

If the lower limit is set to Z m i„=3, for comparison with 
other works, the maximum value, 7.1 fragments, is reached 
at 50 AMeV. The same maximum value is obtained for the 
asymmetric Ni+ Au collisions measured with INDRA (3^ . 
These results can be compared with a series of measure- 
ments of fragment multiplicities published in 32 . Using 
for the present data a similar selection method would in- 
crease the INDRA value by ~10%. These INDRA data, 
combined with those related to Ni+Ni and Kr+Nb of |32~] 
show that the maximum average number of fragments pro- 
duced in central collisions (multifragmentation) is pro- 
portional to the total system mass A sys up to at least 
A sys ^250; there is no indication of a saturation, as in- 
ferred from the Kr+Au data of [33] and from percola- 
tion calculations. The maximum number of fragments for 
A sys ~250 is moreover reached for a cm. available energy 
of ~ 12 AMeV in the INDRA measurements, much lower 
than those quoted in |32|. 

Pictures f), g), h) in the right column of fig. ^show 
the distributions of the sum of the charges of the reaction 
products having atomic numbers Z > 3, 5, and 8, respec- 
tively. The independence of the energy, seen in fig. ^), is 
progressively removed from top to bottom and the pop- 
ulated domain shifts towards lower values of the total 
charge of the reaction products taken into account. For 
Z > 5, for example, the shapes of the distribution are 
very alike: gaussians which are slightly broadening with 
increasing energy, but centred at lower and lower values 
of the sum Abound of the fragment charges. For Z > 8 and 
the largest incident energies (45 and 50 AMeV) shoulders 
are observed at low Abound, which indicates the onset of 
events with a high degree of fragmentation: those events 
have only one fragment with a charge in the range 8-14. 

The evolution of the exit channel average measured 
quantities with the available energy is also synthesized in 
Tabled These average values concern: the total multiplic- 
ity, M tot , the fragment (Z > 5) multiplicity, Mf, the total 
charge emitted in fragments, Abound, and the fragment 
atomic number, Zf. The growth of the average total mul- 
tiplicity, when the incident energy increases, is accounted 
by the growth of the light nuclei (Z < 4) average mul- 
tiplicities, while the average fragment multiplicity barely 
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Fig. 1. Experimental multiplicity (left) and total charge (right) distributions for the selected events from central 25, 32, 39, 
45 and 50 AMeV 129 Xe + nat Sn collisions, with different lower thresholds Zmin for the reaction products taken into account. 
From top to bottom: Z min =1, 3, 5, 8. Note that panel e) thus displays the total detected charge. 



changes. The ratio of the bound charge to the total charge 
of the system goes down with increasing excitation energy. 
The size of the produced fragments is diminishing too, as 
a proof that the collisions are becoming more violent. At 
a given incident energy, however, Abound increases roughly 
linearly as a function of fragment (Z> 5) multiplicity with 
a slope varying from 3 to 4.5 when the energy increases 
from 25 to 50 AMeV. This again indicates that at higher 
energy the system is highly fragmented. 

3.2 Charge distributions 

The charge distributions, normalized to the number of sin- 
gle source events, can be compared all together in fig. 
When the available energy increases, the emission of small 
fragments (Z < 10) increases while that of fragments with 



Z > 15 decreases producing steeper and steeper distribu- 
tions; quite equivalent rates of production for fragments 
with 10 < Z < 15 are observed between 32 and 50 AMeV. 
The domain of the heaviest fragments Z 50, populated 
at the lowest energy, gradually vanishes at higher ener- 
gies. This behaviour is quite remarkable, particularly see- 
ing that the number of fragments is the same at all ener- 
gies (fig. (IJ;). Roughly, one observes the following trend: 
the curves are gradually passing from a regime with two 
slopes, at the lowest incident energy, towards a regime 
with one slope, at the highest incident energy. The av- 
erage measured charge of the fragments can be found in 
Tabled 
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Fig. 2. Experimental differential charge multiplicity distribu- 
tions for the selected events formed in central 25, 32, 39, 45 
and 50 AMeV 129 Xe + nat Sn collisions. 



4 Charge correlation functions and enhanced 
production of nearly equal-sized fragment 
events 

A few years ago a new method called higher order charge 
correlations was proposed to enlighten any extra pro- 
duction of events with specific fragment partitions. The 
high sensitivity of the method makes it particularly ap- 
propriate to look for small numbers of events as those ex- 
pected to have kept a memory of spinodal decomposition 
properties. Thus, such a charge correlation method allows 
to examine model independent signatures that would in- 
dicate a preferred decay into a number of equal-sized frag- 
ments in events from experimental data or from simula- 
tions. In this section and the following only events with 3 
to 6 fragments (Z > 5) at 32, 39, 45 and 50 AMeV will 
be considered. They represent in all cases 90-92 % of the 
selection previously considered. At these four energies, the 
detection efficiency is independent of the fragment charge, 
in the selected event samples. This is not so well verified 
at 25 AMeV. 



4.1 Methods 

The classical two fragment charge correlation method con- 
siders the coincidence yield Y(Z\, Z2) of two fragments of 
atomic numbers Zi t 2, in the events of multiplicity Mf of a 
sample. A background yield Y'(Z\,Z2) is constructed by 
mixing, at random, fragments from different coincidence 
events selected by the same cut on Mf. The two particle 
correlation function is given by the ratio of these yields. 
When searching for enhanced production of events which 
break into equal-sized fragments, the higher order correla- 
tion method appears much more sensitive. All fragments of 



one event with fragment multiplicity Mf = M = nz, 
where nz is the number of fragments with charge Z in 
the partition, are taken into account. By means of the 
normalized first order: 



and second order: 



z 



(1) 



(2) 



moments of the fragment charge distribution in the event 
one may define the higher order charge correlation func- 
tion: 



1 



(3) 



Here, the numerator Y{pz, (%)) is the yield of events 
with given (Z) and az values. Because the measurement of 
the charge belonging to a given event is not subject to sta- 
tistical fluctuations, we use here expression J2J) rather than 
the "nonbiased estimator" of the variance, M _ x ^ z nz{Z— 
(Z)) 2 , as proposed in and used in our previous pa- 
per [TJ\. Note that this choice has no qualitative influ- 
ence on the forthcoming conclusions. The denominator 
Y'(az } (Z)), which represents the uncorrelated yield of 
pseudo-events, was built in 34, as for classical correla- 
tion methods, by taking fragments at random in different 
events of the selected sample of a certain fragment multi- 
plicity. This Monte-Carlo generation of the denominator 
Y' (az, (Z)) can be replaced by a fast algebraic calculation 
which is equivalent to the sampling of an infinite number 
of pseudo-events ^Sj. Its contribution to the statistical 
error of the correlation function is thus eliminated. How- 
ever, owing to the way the denominator was constructed, 
only the fragment charge distribution dM/dZ of the par- 
ent sample is reproduced but the constraints imposed by 
charge conservation are not taken into account. This has, 
in particular, a strong effect on the charge bound in frag- 
ments dAf/dZbound distribution. This fact makes the de- 
nominator yield distributions as a function of (Z) wider 
and flatter than those of the numerator |2H]. Consequently, 
even in the absence of a physical correlation signals, the 
ratio © is not a constant equal to one. The correlations 
induced by the finite size of the system (charge conser- 
vation) distorts the amplitude, or even may cancel other 
less trivial correlations. Therefore, a new method for the 
evaluation of the denominator ^B] , based on the "intrinsic 
probability" of emission of a given charge, was adopted. 
It minimizes these effects and replicates all features of the 
partitions of the numerator, except those (of interest) due 
to other reasons than charge conservation. 

The goal of the method is to take into account in a 
combinatorial way the trivial correlations due to charge 
conservation. If there is no correlation between the charges, 
each charge can be fully described by an emission probabil- 
ity (referred to as intrinsic probability and noted mtl 'Pz). 
This fact was largely demonstrated for the bulk of multi- 
fragment emission around the Fermi energy |35j . Without 
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charge conservation constraint, the intrinsic probabilities 
would be equal to the emission probability. Charge con- 
servation makes the emission probability of high charge 
fragments smaller than their intrinsic probability, whereas 
it is the contrary for small charges. The probability to ob- 
serve a given partition (n : (n\, . . . , riz tot )), at a given 
total multiplicity, M to t — m — 53 z n Zi is obtained by the 
multinomial formula. If the total charge is fixed (Z to t = 
Z nz), the partition probabilities are given by: 



intr r>nz 

P(n\m) =am\ R-^f- 6z tot ,E z z nz , (4) 



where a is the normalization constant (so that ^ n P(n\m) 
= 1) and S is the Kronecker symbol. For the reason previ- 
ously mentioned, no Z-dependent efficiency term is needed 
in this formula. These partition probabilities contain all 
information relative to the charges and their correlations. 
For example, the denominator Y' (az, (Z)) is obtained by 
summing the probabilities of the partitions with given 
mean charge and standard deviation. Of course, the intrin- 
sic probabilities are not direct experimental observables, 
they have to be evaluated by inversion of Eq. Q . However, 
this inversion is possible only if the physical correlations 
(not due to charge conservation) are weak. If the data sam- 
ple contains only trivial correlations, then the higher order 
charge correlation function (or any other correlation func- 
tion) is everywhere equal to 1. Local positive (negative) 
physical correlations will appear as peaks (holes). In the 
case of experimental events, it has to be noted that a set 
of intrinsic probabilities exists only if the selected sample 
comprises essentially single source events having reached 
thermal equilibrium (the source has to be described by 
a unique set of intrinsic probabilities). Otherwise the re- 
sults would be equal to the convolution of the partition 
probabilities corresponding to the different source sizes, 
each of them being described by a different set of intrin- 
sic probabilities. Hence, the convergence of the inversion 
procedure (Eq. 0} is a strong indication that these two 
conditions are fulfilled by the sample under study |18| . 
The data samples correspond to multifragmenting single 
sources, supplemented by pre-equilibrium and secondary 
decay particles and light fragments. Hence, only fragments 
with charge greater than or equal to 5 have been retained 
for the analysis of the multifragmenting source. This fea- 
ture has to be taken into account in the calculation of the 
denominator. The fragment partition (N : (72,5, . . . , nz tat )) 
probability can be calculated as the sum over the complete 
partitions n which include N. Noting M the fragment mul- 
tiplicity, it can be straightforwardly shown from Eq. 
that the fragment partition probability reads: 



P(N) =/(M,Z bound )P'(N) , (5) 



with 

f(M, Z hound ) =a £ (m) ( m " M ) ! 

n:(ni,...,n 4 ) 

4 intrpnz ( 6 ) 

X 11 n7 \ 6 Z tat -Z bouad ,T.% =1 Zn z 
Z=\ 

P'(N) = Ml U Z z= 5 , 

where (^) are the binomial coefficients. The value of 
Z to t was fixed at the total charge of the system, Z pro jectiie 
+ -^target; it has been noticed that, whereas the total 
charge conservation has to be explicitely included in the 
calculation, the results are only modified when Z to t is 
lower than 85. This new method to build the denomina- 
tor will be denoted as the Intrinsic Probability Method 
(IPM) in what follows. However, the explicit calculation 
of the intrinsic probabilities may not be the only method 
for building a denominator including only the correlations 
induced by charge conservation (another procedure is pro- 
posed in see Appendix for a comparison with the IPM 
method). In any case, by definition, all valid methods lead 
to a partition probability described by Eq. QJ. 

A comparison of the results obtained in the framework 
of the higher order correlation method, with the analytical 
denominator, and with IPM, will be drawn in the follow- 
ing for an event sample resulting from a simulation of 32 
AMeV 129 Xe + 119 Sn central collisions. 



4.2 Stochastic mean-field and spinodal decomposition 

Dynamical stochastic mean-field simulations have been 
proposed for a long time to describe processes involv- 
ing instabilities like those leading to spinodal decomposi- 
tion 37 38 39 . In this approach, spinodal decomposition 
of hot and dilute finite nuclear systems can be mimicked 
through the Brownian One-Body (BOB) dynamics |4(JI 
I41IIT] . which consists in employing a Brownian force in 
the kinetic equations. Simulations have been performed 
for head-on 129 Xe on 119 Sn collisions at 32 AMeV. The 
ingredients of simulations can be found in |2j as well as 
a detailed comparison between filtered simulated events 
(to account for the experimental device) and experimen- 
tal data. A good agreement between both is revealed. 

To refine the comparison higher-order charge correla- 
tions have been calculated for the simulated events, keep- 
ing the compact presentation proposed in |17j : charge cor- 
relation functions are built for all events, whatever their 
multiplicity, by summing the correlated yields for all M 
and by replacing the variable (Z) by Abound = M x (Z) = 

Znz- Uncorrelated events are constructed and weighted 
in proportion to real events of each multiplicity. This pre- 
sentation is based on the experimental observation that 
the peaks observed independently for each fragment mul- 
tiplicity correspond to the same Abound region rT7|. The 
variance bin was chosen equal to one charge unit. We re- 
call that in the considered domain of excitation energy, 
around 3 MeV per nucleon |3()ll7j . secondary evaporation 
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Fig. 3. Correlation functions for events with Mf = 3 to 6, simulated with the Brownian One Body model for 32 AMeV 1 Xe 
+ nat Sn collisions, a) with an analytical denominator provided by pseudo-events; b) with a denominator calculated with the 
IPM. The orientations of a) and b) are different for a better visualisation of the landscapes 



leads to fragments one charge unit smaller, on average, 
than the primary Z rj 10 — 20 ones, with a standard de- 
viation around one |29| . If a weak enhanced production of 
exactly equal-sized fragments exists, peaks are expected 
to appear in the interval az = — 1, because of secondary 
evaporation. This interval in az is hence the minimum 
value which must be chosen to look for nearly equal-sized 
fragments. Any (unknown) intrinsic spread in the frag- 
ment size coming from the break-up process itself may 
enlarge the az interval of interest. This will be discussed 
at the end of the section, for the moment we will consider 
only events with az < 1, which corresponds to differences 
of at most two units between the fragment atomic num- 
bers in one event. 

Fig. shows the correlation function calculated us- 
ing the analytical denominator (a) or the denominator 
given by the IPM (b). Both functions are drawn versus 
the variables Abound = M x (Z) and az- In fig. Ek, the 
equal-sized fragment correlations in the first bin are su- 
perimposed over trivial correlations due to the finite size 
of the system. For this reason, the ratio is gener- 
ally different from one and smoothly varies with the vari- 
ables Abound and az- For each bin in Zbound (fixed at 
6 atomic number units), an exponential evolution of the 
correlation function is observed from az — 7 — 8 down to 
az = 2 — 3. This exponential evolution is thus taken as a 
"background" empirically extrapolated down to the first 
bin az = — 1 The amplitude of the correlation function 
in the domain Abound = 36 — 60 is well above the back- 
ground, with a confidence level higher than 90%, proving 
thus a statistically significant enhancement of equal-sized 
fragment partitions. Of the 1% of events having az < 1, 
(0.13±0.02)% (called extra-events from now on) are in ex- 
cess of the background. In fig. |2h>, as expected, all cor- 



relations due to the charge conservation are suppressed 
and the correlation function is equal to 1 (within statis- 
tical fluctuations) wherever no additional correlation is 
present. Again one observes peaks for az < 1. The per- 
centage of extra-events is 0.36±0.03%, higher than the 
one obtained with the previous method. Moreover, with 
this method, peaks also appear at the maximum values 
of az for a given Zbound- They correspond to events com- 
posed of one big (a heavy residue) and several lighter frag- 
ments (sequentially emitted from the big one) . In that case 
fusion-multifragmention does not occur and the peaks re- 
veal the small proportion (0.15%) of events which undergo 
the fusion-evaporation process. The IPM approach may 
thus reveal other correlations not seen with the previous 
one. Therefore, we shall use this approach in the following 
analysis of the experimental data. 

To conclude this part we can say that, although all 
events in the simulation arise from spinodal decompo- 
sition, only a very small fraction of the final partitions 
have nearly equal-sized fragments. Different effects: beat- 
ing of different modes, coalescence of nascent fragments, 
secondary decay of the excited fragments and, above all, 
finite size effects arc responsible for this fact [151116] . The 
signature of spinodal decomposition can only reveal itself 
as a "fossil" signal. 

4.3 Experimental results 

We shall present now higher-order charge correlations for 
the selected experimental events. This will be done for 
four incident energies (32, 39, 45 and 50 AMeV) in the 
framework of the IPM for the denominator. The first step 
consists in determining the intrinsic probabilities of frag- 
ments for each multiplicity and at each incident energy. 
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Fig. 4. Experimental differential charge multiplicity distribu- 
tions (circles) for the single source formed in central 39 AMeV 
129 Xe on nat Sn collisions. Parts a, b, c and d refer respectively 
to fragment multiplicities 3, 4, 5, 6. The Z distributions for the 
first (squares), second (diamonds) and third (triangles) heav- 
iest fragments are presented too. The lines correspond to the 
results obtained with IPM. The dashed lines display the intrin- 
sic probabilities. 



These probabilities are obtained by a recursive procedure 
of minimization. The minimization criterion is the normal- 
ized x 2 between the experimental fragment probabilities 
and the fragment partition probabilities given by JSJ- The 
calculated x 2 were always lower than one: the lower is the 
incident energy and the larger is the multiplicity, the lower 
is the x 2 value. 

Charge distributions experimentally observed for the 
different fragment multiplicities are shown as an example 
in fig. they correspond to the 39 AMeV incident energy. 
Dashed lines refer to the intrinsic probabilities calculated 
with IPM and the corresponding charge distributions are 
the full lines. The experimental charge distributions are 
faithfully described. Note that the ratio of the experimen- 
tal to the calculated curves does not reveal any anomalous 
enhancement in a preferential domain of fragment atomic 
number. The charge distributions, summed over fragment 
multiplicities 3 to 6, are displayed in fig. [H] for the four 
incident energies. For each incident energy the intrinsic 
probabilities have been calculated independently for the 
different fragment multiplicities (see fig. 0} and weighted 
in proportion to real events of each fragment multiplicity. 
The slight differences for the intrinsic probabilities corre- 
sponding to the different multiplicities reveal small differ- 
ences in the average excitation energies of the multifrag- 
menting sources (on the other hand, it has been possible 
to fit all events from the BOB simulation, irrespective of 



Fig. 5. Experimental differential charge multiplicity distribu- 
tions (circles) for the single sources formed in central 32, 39, 45 
and 50 AMeV 129 Xe on nat Sn collisions (a, b, c, d). The Z dis- 
tributions for the first (squares), second (diamonds) and third 
(triangles) heaviest fragments are presented too. The lines cor- 
respond to the results obtained with IPM. 



their multiplicity, using the same intrinsic probability dis- 
tribution due to the fact that the events correspond to 
the same initial conditions). We note again the excellent 
agreement between calculations and data at all incident 
energies. The contributions to the Z distribution of the 
three heaviest fragments of each partition are as well de- 
scribed as in fig. 01 and the charges bound in fragments 
(not shown) are also perfectly reproduced. 

Fig. illustrates the higher-order correlation functions 
measured for the different fragment multiplicities. It con- 
cerns single sources selected at 39 AMeV incident energy. 
To make the effects more visible, peaks with confidence 
level lower than 80% were flattened out. We observe sig- 
nificant peaks in the bin a — — 1 for each fragment 
multiplicity. For M=6, peaks are essentially located in the 
bin oz = 1 — 2. As observed in simulations, peaks corre- 
sponding to events composed of a heavy residue and light 
fragments (oz in the region 5-10 associated with low (Z)) 
are also visible. 

Let us now present the results for the different inci- 
dent energies, summed over multiplicities 3-6. We recall 
again that IPM denominators of correlation functions are 
weighted in proportion to real events of each multiplicity. 
In fig. the measured functions are displayed. If we ex- 
clude some peaks with low confidence levels (as done in 
fig. 01 correlation functions are equal or close to one except 
at low oz values. A summary of average-charge domains 
contributing to the correlation peaks in the first bin in oz 
is given in Table as a function of fragment multiplicity. 
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Fig. 6. Experimental higher-order charge correlations for selected events formed in central 39 AMeV 129 Xe on nat Sn collisions, 
for fragment multiplicities 3 to 6. The maximum value of the scale of the correlation function is limited to 3 on the picture. 



Table 2. Characteristics of events with az < 1 for the dif- 
ferent incident energies. For each fragment multiplicity M, the 
range of (Z) contributing to the correlation peaks are indi- 
cated (for bold (Z) ranges, see text). The first line refers to 
BOB simulations (see subsect. I4.2H . 



M 

E (AMeV) 


3 






4 




5 


6 


32* 


12 - 


20 


9 


- 17 


8 


- 13 




32 


13 - 


21 


11 


- 16 








39 


15 - 


20 


8 


- 15 


7 


- 11 


6 - 8 


45 


17 - 


18 


5 


- 14 


8 


- 11 


6 - 8 


50 














7 - 9 



All multiplicities, associated to the largest (Z) ranges (20- 
6 and 18-6) contribute to the peaks at 39 and 45 AMeV. At 
32 AMeV the peaks only come from the low multiplicites 
(M = 3-4), with a smaller (Z) range : 21-11; at 50 AMeV, 
the situation is completely different and we only observe 
a contribution from M = 6 with (Z) = 7-9. Except in this 
latter case, the (Z) domains for which correlation peaks 
are present are similar but slightly shifted towards lower 
values with increasing energy. 



To be quantitative, correlation functions for the first 
bin in az are displayed in fig. with their statistical er- 
rors. Full points (open points) correspond to a confidence 
level higher (lower) than 90%. Two percentages of events 
relatively to the single source events (Mf = 3 — 6) are pre- 
sented in columns 2 and 3 of Tabled they refer to the total 
number of events and of extra-events (taking into account 
correlations and anti-correlations) in the bin. The present 
analysis fully confirms at 32 AMeV incident energy the 
previous one ^ZJ and the extra-percentage of events with 
nearly equal-sized fragments is maximum at 39 AMeV. 

A closer examination of fig. Eh and[H]reveals an islet of 
peaks, with a high confidence level, in the second and the 
third bins in az- They arc located in the upper region of 
the M x (Z) ranges indicated in table El They correspond 
to events with a broader spread of charges. Observing such 
events in the simulation may indicate that the intrinsic 
spread (which is unknown) due to spinodal decomposition 
is larger than that coming from secondary evaporation, 
and hence that these peaks also sign the original process. 
The percentages of events with az < 3 are also reported 
in table El The conclusions are the same as above: while 
more events have small values of az when the incident 
energy increases, the percentage of extra-events shows a 
maximum at 39 AMeV, and vanishes at 50 AMeV. 
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Fig. 7. Experimental higher-order charge correlations for selected events formed in central 129 Xe and nat Sn collisions. Events 
with M/=3 to 6 are mixed. 



Table 3. Numbers and percentages of events for az < 1 and oz < 3 and percentages of extra-events. The percentages 
are expressed with respect to the total number of selected events with fragment multiplicities 3 to 6. The first line refers to 
simulations (see subsect. l4.2H . 

oz < 1 o~z < 3 



E 


events 


extra-events 


events 


extra-events 


(AMeV) 




(%) 


(%) 




(%) 


(%) 


32* 


353 


1.0 


0.36 ± 0.03 


4873 


14.0 


1.30 ± 0.06 


32 


83 


0.26 


0.13 ± 0.02 


1746 


5.5 


0.79 ± 0.05 


39 


151 


0.58 


0.25 ± 0.03 


4275 


16.3 


1.29 ± 0.07 


45 


317 


1.32 


0.21 ± 0.03 


7328 


30.6 


0.79 ± 0.05 


50 


762 


2.78 


0.08 ± 0.02 


12306 


45 






5 Discussion 

The first hint of a bulk effect for producing fragments in 
central collisions between heavy nuclei at a given mod- 
erate excitation energy (e* 7 AMeV) was the measure- 
ment of identical fragment charge distributions for two 
different system sizes (32 AMeV 129 Xe + ™ at Sn and 36 
AMeV 155 Gd + nat U) g2 However this feature could 
as well be interpreted from a statistical point of view, 
namely the dominance of phase space |35|43) . Indeed both 
dynamical and statistical approaches were able to repro- 
duce the experimental observation. With the same dy- 
namical approach as that used in this paper, the experi- 
mental charge and multiplicities distributions were repro- 
duced while the average fragment kinetic energies were 



accounted for within 20 % [7] . These same properties were 
also well accounted for with the statistical model SMM |431 
44 45 , 46 . From this agreement between data and the two 
models we learnt that the dynamics involved is sufficiently 
chaotic to finally explore enough of the phase space and 
describe fragment production through a statistical approach. 
To go further more constrained observables were needed. 
Such was the goal of the present studies. 

We have just seen that, at 32 AMeV incident energy, 
experimental correlation functions are similar to those ob- 
tained with events from dynamical simulations, BOB. Above 
all, in both of them, there is abnormal production of nearly 
equal-sized fragments pointed out by the peaks in the first 
oz bin(s). Supported by this simulation we thus attributed 
the greatest part of fusion-multifragmentation events to 
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Fig. 8. Same as fig.[7|but peaks (or holes) with a confidence level lower than 90% (1.65 a) have been flattened out. 



spinodal decomposition J7| . These peaks hence appear as 
fossil fingerprints of the break-up process. 

Correlation peaks are also observed for experimental 
data at higher incident energies: the excitation function is 
displayed in fig. 1101 Information on the associated thermal 
excitation energies (and extra radial collective energy) in- 
volved over the incident energy domain studied can be pro- 
vided by the SMM model which well describes static and 
dynamics observables of fragments. Starting from a freeze- 
out volume fixed at three times the normal volume, the 
thermal excitation energies of the dilute and homogeneous 
system, extracted from SMM, vary from 5.0 to 7.0 AMeV 
and the added radial expansion energy remains low: from 
0.5 to 2.2 AMeV |33I37| . A rise and fall of the percentage 
of "fossil partitions" from spinodal decomposition is mea- 
sured. Fig. 1101 reveals some difference between the exper- 
imental (full symbols) and simulated events (open sym- 
bols): the experimental percentages of extra events are 
closer to the simulated ones in fig. 1 10b than in fig. 110b . 
This means that the charge distributions inside an event 
are slightly narrower in the simulation than in the exper- 
iment either because of the primary intrinsic spread, or 
because the width due to evaporation is underestimated. 
For the considered system, incident energies around 35- 
40 AMeV appear as the most favourable to induce spin- 
odal decomposition; it corresponds to about 5.5-6 AMeV 
thermal excitation energy associated to a very gentle ex- 
pansion energy around 0.5-1 AMeV. The qualitative ex- 
planation for those numbers can be well understood in 



terms of a necessary compromise between two times. On 
one hand the fused systems have to stay in the spinodal 
region w 100-150 fm/c 16,48,49 , to allow an important 
amplification of the initial fluctuations and thus permit 
spinodal decomposition; this requires a not too high inci- 
dent energy, high enough however for multifragmentation 
to occur. On the other hand, for a finite system, Coulomb 
interaction and collective expansion push the "primitive" 
fragments apart and reduce the time of their mutual inter- 
ation, which is efficient to keep a memory of "primitive" 
size properties. Above 45 AMeV incident energy part of 
trajectories followed by the system in the temperature- 
density plane may not sufficiently penetrate the spinodal 
region. Finally, considering secondary decay effects, they 
should be essentially similar between 39 and 50 AMeV be- 
cause the primary fragment excitation energy per nucleon 
was found constant over this incident energy range |27| . 
Particularly the excitation function shown in fig.llOlshould 
not be affected by such effects. 

Let us come now to the size of fragments associated 
to "fossil partitions" . From the theoretical point of view, 
spinodal instabilities have been mainly studied within semi- 
classical (as BOB) or hydrodynamical |49II50II51| frame- 
works. Precise information on the most unstable collective 
modes, which is needed for discussing fragment size pro- 
duced by spinodal decomposition, can only be obtained 
when considering quantal effects. From the limited num- 
ber of studies investigating such effects [T4ll52ll53| a few 
trends emerge: i) an increase of the number of unstable 
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Fig. 10. Abnormal production of events with nearly equal-sized fragments (a: az < 1 and b: az < 3) as a function of thermal 
excitation energy (deduced from SMM): full points. The incident and radial energy scales are also indicated. The open point 
refers to the result from BOB simulations; the average thermal excitation energy is used. Vertical bars correspond to statistical 
errors and horizontal bars refer to estimated uncertainties on the backtraced quantity, e*. 
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Fig. 9. Higher-order charge correlations at the different in- 
cident energies: quantitative results for events with az < 1. 
Errors are statistical assuming independent measurements and 
horizontal bars correspond to the M x (Z) bins. 



collective modes with the size of the system (up to mul- 
tipolarity L—5-6 for A=140), ii) a decrease of the L max 
value when the temperature increases and iii) a dominance 
of the octupole mode (L=3) which appears as the most 
unstable. Charge correlations were only studied for events 
with 3 to 6 fragments, which precludes any information on 



the quadrupole mode (L=2). If we except results at the 
highest incident energy, which strongly differ from oth- 
ers, we do observe the dominance of events with three 
fragments, associated to the multipolarity L=3. Indeed 
the Z domains marked in bold in table [5] correspond to 
the largest proportion of events which populate the first 
a z bin, which vary from 70% at 32 AMeV to 40% at 45 
AMeV. Then the number of events with higher multiplic- 
ities, corresponding to larger L values, progressively de- 
creases, as expected. 

At 50 AMeV results observed for the fragment sizes 
(az < 1) are puzzling. Moreover the percentage of extra- 
events drops to zero when events with az < 3 are consid- 
ered. Further simulations and theoretical works are needed 
to progress in the interpretation of these data. 



6 Conclusions 

In conclusion, we have investigated charge correlation func- 
tions for compact and heavy fused systems which undergo 
multifragmcntation, as a function of the incident energy, 
from 32 to 50 AMeV. At the lowest energy, we have con- 
firmed an enhanced production of events with equal-sized 
fragments. Supported by theoretical simulations we have 
interpreted this enhancement as a signature of spinodal 
instabilities as the origin of multifragmentation of those 
systems in the Fermi energy domain. This fossil signal 
culminates for incident energy around 35-40 AMeV, which 
corresponds to the formation of a hot and dilute system 
at 0.3-0.4 the normal density and temperature around 4- 
5 MeV |5~4lT7| . Spinodal decomposition describes the dy- 
namics of a first order phase transition, and the present 
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observations support the existence of such a transition for 
hot finite nuclear matter [3"ll45ll9"] . 

The experimental results were obtained by the INDRA collab- 
oration. 



Appendix 
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Fig. 11. Deviations from 1 of the correlation functions divided 
by the statistical errors in abcissa, for the different values of 
az- Correlation functions calculated by the IPM method, and 
shown in fig0 



Very recently a new method for building charge cor- 
relation functions was proposed, of implementation easier 
than the IPM method: the denominator is built by mix- 
ing events through random exchanges of two fragments 
between two events under the constraint that the sum 
of the two exchanged fragments is conserved, besides the 
conservation of Abound EH]- To test the sensitivity ot this 
method as compared to the IPM method, the present ex- 
perimental data were also analyzed using the exchange 
method, by building correlation functions 1 + R(az, {Z}) 
with a z being calculated from eq.|3 The results of the two 
methods are compared in fig. II II and fig. 1121 displaying for 
each bin of the plane (az, M x(Z)), the deviation from 
1, R(az, (Z)), of the correlation function, normalized to 
its statistical error bar, <Ji+R( az .(z)), calculated from the 
numerator only. 

The greater sensitivity of the IPM method is clearly 
seen on fig. ^2 the values of the correlation function are 
closer to one (R=0) except at low az where we observe 
correlations with a significant confidence level. Conversely, 
the exchange method, fig. El leads to a large dispersion of 
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Fig. 12. Deviations from 1 of the correlation functions divided 
by the statistical errors in abcissa, for the different values of 
az- Correlation functions calculated by the exchange mixing 
event method. 



the values of R(az, (•Z))/o'i+h( - z ,(z))! ~1.6 times broader 
than with the IPM method at 32, 39 and 45 AMeV. At 
50 AMeV both methods lead to similar dispersions. 



References 

1. L. G. Moretto and G. J. Wozniak, Ann. Rev. of Nuclear 
and Particle Science 43, 379 (1993). 

2. D. H. E. Gross, Rep. Prog. Phys. 53, 605 (1990). 

3. M. D'Agostino, F. Gulminelli, P. Chomaz, M. Bruno, 
F. Cannata et al, Phys. Lett. B 473, 219 (2000). 

4. J. A. Hauger et al. (EOS collaboration), Phys. Rev. C 62, 
024616 (2000). 

5. J. B. Elliott et al. (EOS collaboration), Phys. Rev. C 62, 
064603 (2000). 

6. H. Xu et al, Phys. Rev. Lett. 85, 716 (2000). 

7. J. D. Frankland, B. Borderie, M. Colonna, M. F. Rivet, 
C. O. Bacri et al. (INDRA collaboration), Nucl. Phys. A 
689, 940 (2001). 

8. R. Botet, M. Ploszajczak, A. Chbihi, B. Borderie, D. Du- 
rand et al, Phys. Rev. Lett. 86, 3514 (2001). 

9. M. D'Agostino, R. Bougault, F. Gulminelli, M. Bruno, 
F. Cannata et al., Nucl. Phys. A 699, 795 (2002). 

10. A. H. Raduta and A. R. Raduta, Phys. Rev. C 65, 054610 
(2002). 

11. M. Kleine Berkenbrusch, W. Bauer, K. Dillman, S. Pratt, 
L. Beaulieu et al, Phys. Rev. Lett. 88, 022701 (2002). 

12. J. B. Elliott, L. G. Moretto, L. Phair, G. J. Wozniak, 
L. Beaulieu et al, Phys. Rev. Lett. 88, 042701 (2002). 

13. G. F. Bertsch and P. J. Siemens, Phys. Lett. B 126, 9 
(1983). 

14. S. Ayik, M. Colonna and P. Chomaz, Phys. Lett. B 353, 
417 (1995). 



14 



G. Tabacaru et al.: Fragment Charge Correlations and Spinodal Decomposition in Finite Nuclear Systems 



15. B. Jacquot, S. Ayik, P. Chomaz and M. Colonna, Phys. 
Lett. B 383, 247 (1996). 

16. M. Colonna, P. Chomaz and A. Guarnera, Nucl. Phys. A 
613, 165 (1997). 

17. B. Borderie, G. Tabacaru, P. Chomaz, M. Colonna, 
A. Guarnera et al. (INDRA collaboration), Phys. Rev. 
Lett. 86, 3252 (2001). 

18. P. Desesquelles, Phys. Rev. C 65, 034604 (2002). 

19. J. D. Frankland, C. O. Bacri, B. Borderie, M. F. Rivet, 
M. Squalli et al. (INDRA collaboration), Nucl. Phys. A 
689, 905 (2001). 

20. J. Pouthas, B. Borderie, R. Dayras, E. Plagnol, M. F. Rivet 
et al, Nucl. Instr. and Meth. in Phys. Res. A 357, 418 
(1995). 

21. J. Pouthas, A. Bertaut, B. Borderie, P. Bourgault, B. Ca- 
han et al., Nucl. Instr. and Meth. in Phys. Res. A 369, 
222 (1996). 

22. M. P&rlog, B. Borderie, M. F. Rivet, G. Tabacaru, 
A. Chbihi et al. (INDRA collaboration), Nucl. Instr. and 
Meth. in Phys. Res. A 482, 674 (2002). 

23. M. Paxlog, B. Borderie, M. F. Rivet, G. Tabacaru, 
A. Chbihi et al. (INDRA collaboration), Nucl. Instr. and 
Meth. in Phys. Res. A 482, 693 (2002). 

24. D. Cugnon and D. L'Hote, Nucl. Phys. A 397, 519 (1983). 

25. J. F. Lecolley, L. Stugge, M. Aboufirassi, A. Badala, B. Bil- 
wes et al., Phys. Lett. B 325, 317 (1994). 

26. N. Marie, R. Laforest, R. Bougault, J. P. Wieleczko, D. Du- 
rand et al. (INDRA collaboration), Phys. Lett. B 391, 15 

(1997) . 

27. S. Hudan, A. Chbihi, J. D. Frankland, A. Mignon, A. Botv- 
ina et al. (INDRA collaboration), Phys. Rev. C, in press 
(2003). 

28. Anh-Dung Nguyen, these de doctorat, Universite de Caen 

(1998) , LPCC T 98-02. 

29. G. Tabacaru, these de doctorat, Universite Paris-XI Orsay 
(2000), IPNO T 00-13. 

30. N. Marie, A. Chbihi, J. Natowitz, A. Le Fevre, S. Sa- 
lou et al. (INDRA collaboration), Phys. Rev. C 58, 256 
(1998). 

31. N. Bellaize, O. Lopez, J. P. Wieleczko, D. Cussol, G. Auger 
et al. (INDRA collaboration), Nucl. Phys. A 709, 367 
(2002). 

32. D. Sisan, W. Bauer, O. Bjarki, D. J. Magestro, A. Nadasen 
et al., Phys. Rev. C 63, 027602 (2001). 

33. C. Williams, W. G. Lynch, C. Schwarz, M. B. Tsang, W. C. 
Hsi et al., Phys. Rev. C 55, 2132 (1997). 

34. L. G. Moretto, T. Rubehn, L. Phair, N. Colonna, G. J. 
Wozniak et al., Phys. Rev. Lett. 77, 2634 (1996). 

35. L. G. Moretto, R. Ghetti, L. Phair, K. Tso and G. J. Woz- 
niak, Phys. Rep 287, 249 (1997). 

36. J. L. Charvet, R. Dayras, D. Durand, O. Lopez, L. Nalpas 
et al., to be published nucl-ex/0212012 

37. J. Randrup and B. Remaud, Nucl. Phys. A 514, 339 
(1990). 

38. P. Chomaz, G. F. Burgio and J. Randrup, Phys. Lett. B 
254, 340 (1991). 

39. G. F. Burgio, P. Chomaz and J. Randrup, Nucl. Phys. A 
529, 157 (1991). 

40. P. Chomaz, M. Colonna, A. Guarnera and J. Randrup, 
Phys. Rev. Lett. 73, 3512 (1994). 

41. A. Guarnera, P. Chomaz, M. Colonna and J. Randrup, 
Phys. Lett. B 403, 191 (1997). 



42. M. F. Rivet, C. O. Bacri, B. Borderie, J. D. Frankland, 
M. Assenard et al. (INDRA collaboration), Phys. Lett. B 
430, 217 (1998). 

43. J. Bondorf, A. S. Botvina, A. S. Iljinov, I. N. Mishustin 
and K. Sneppen, Phys. Rep. 257, 133 (1995). 

44. S. Salou, these de doctorat, Universite de Caen (1997), 
GANIL T 97 06. 

45. N. Le Neindre (INDRA collaboration), G. Bonsignori, 
M. Bruno, A. Ventura and D. Vretenar, eds., Proc. of the 
International Conference on Structure of the Nucleus at 
the Dawn of the Century, Bologna, Italy (World Scientific) 
(2001), Nucleus-Nucleus collisions, 221. 

46. G. Tabacaru (INDRA collaboration), G. Bonsignori, 
M. Bruno, A. Ventura and D. Vretenar, eds., Proc. of the 
International Conference on Structure of the Nucleus at 
the Dawn of the Century, Bologna, Italy (World Scientific) 
(2001), Nucleus-Nucleus collisions, 321. 

47. N. Le Neindre, these de doctorat, Universite de Caen 
(1999), LPCC T 99 02. 

48. D. Idier, M. Farine, B. Remaud and F. Sebille, Ann. Phys. 
Fr. 19, 159 (1994). 

49. W. Norenberg, G. Papp and P. Rozmej, Eur. Phys. J. A 
14, 43 (2002). 

50. W. Norenberg, G. Papp and P. Rozmej, Eur. Phys. J. A 
9, 327 (2000). 

51. H. Heiselberg, C. J. Pethick and D. G. Ravenhall, Phys. 
Rev. Lett. 61, 818 (1988). 

52. P. Chomaz, D. Lacroix, B. Jacquot, M. Colonna and 
S. Ayik, H. Feldmeier, J. Knoll, W. Norenberg and 
J. Wambach, eds., Proc. XXVII Int. Workshop on Gross 
properties of nuclei and Nuclear excitations, Hirschegg, 
Austria (1999), 312. 

53. M. Colonna, P. Chomaz and S. Ayik, Phys. Rev. Lett. 88, 
122701 (2002). 

54. A. Guarnera, these de doctorat, Universite de Caen (1996), 
GANIL T 96 01. 



